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Light meson correlation functions near the deconfining transition on 
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We study the light hadron correlators near the deconfining transition by extracting the spectral function on 
quenched anisotropic lattices. We adopt the method successfully applied to the charmonium systems: the use 
of the smeared operators and the analysis method which is composed of the maximum entropy method and the 
fitting assuming several forms. The numerical simulations are performed on lattices of aj^ ~ 2 GeV and 
as/at — 4, with the clover quark around the strange quark mass. Towards applications to the plasma phase, we 
check the reliability of our procedure by examining how the results for the correlators below Tc {Nt — 32) changes 
under variation of input parameters such as the smearing function. 



1. Introduction 

The theoretical expectation for properties of 
quarks and gluons in the deconfinement phase is 
strongly required for an interpretation of exper- 
imental results in the heavy ion collision experi- 
ments. Contrary to a naive quark gluon plasma 
(QGP) picture such as the gas of weakly inter- 
acting quarks and gluons, several model calcula- 
tions inspired by QCD indicate nonperturbative 
features of QGP. For example, an NJL model cal- 
culation suggests that the hadronic excitations 
survive as soft modes even above the critical tem- 
perature pp. 

The lattice QCD can in principle solve the 
problem in a model independent way, since the 
dynamical properties of excitation modes can be 
extracted from the Matsubara Green functions 
[2]. Although in practice such studies are quite 
nontrivial even at qualitative level, recent techni- 
cal development such as anisotropic lattices and 
computational progress have enabled us to study 
dynamical properties of hadrons at finite tem- 
perature |3I4| . In addition, several procedures 
such as the maximum entropy method (MEM) 
[5] have been developed to analyze the spectral 
function with which we can directly probe the ex- 
citation modes of the plasma phase. These proce- 
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dures have been applied to the charmonium sys- 
tems and suggest that hadronic excitations sur- 
vive even above the critical temperature |6l7l8j . 

In this work we apply the method which was 
adopted in Ref. [HI to systems with light quarks. 
In this report we focus on a check of reliability 
for applied analysis methods at finite tempera- 
ture, as a preparation for future applications to 
the deconfinement phase. 



2. Our approach 

Here we briefly explain our analysis procedure 
and its criteria for a reliability, which were in de- 
tail described in Ref. [H]. 

In order to extract reliable information on the 
spectral function, we use MEM and fit method 
in a complementary manner. After a rough esti- 
mation of spectral function by MEM, we eval- 
uate properties of the mode such as mass and 
width more quantitatively by the fit assuming 
several fitting forms based on MEM results. As 
more sophisticated form of the fit we adopt the 
constrained curve fitting (CCF) |S] whose prior 
knowledge is also estimated with the MEM re- 
sults. 

To verify the reliability of our analysis, we 
require the following criteria for the extraction 
methods: (1) The stability of the spectral func- 
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tion against variations of input parameters or 
model functions. (2) The stability of the result 
for the correlators at T = under restriction of 
the degrees of freedom. For the latter, since the 
T = correlators are not ready at present, we 
instead analyze those below Tc which should pro- 
duce similar results as at T = 0. 

In the previous work for the charmonium, in or- 
der to satisfy the criteria we employ the smearing 
operators which enhance the low frequency part 
of the meson correlators. However the smeared 
operators may produce artificial peaks in the 
spectral function ^H). In order to distinguish an 
artificial peak from the genuine physical one we 
apply several smearing functions and examine the 
stability of the results. 
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Figure 1 . tmax dependence of the spectral func- 
tion for the point correlator in PS channel. 



3. Numerical results 

We use quenched lattices of the sizes 20'^ x Nt, 
where Nt — 160 and 32 which roughly corre- 
spond to r ~ and 0.9Tc, respectively. The 
zero temperature lattice, and the setup of lat- 
tice parameters are the same as those used in 
Ref. JI]. The gauge configurations are gener- 
ated with the Wilson plaquette action at /? = 6.10 
with the renormalized anisotropy ^ = a^/at = 4. 
The spatial cutoff set by the hadronic radius vq 
is = 2.030(13) GeV. Nt = 28 is close to the 
phase transition. The quark action is the 0{a) 
improved Wilson action 11 , with the hopping 
parameter roughly corresponding to the strange 
quark mass. We have 500 configurations at T = 
and 1000 configurations at T > 0. 

We measure the point correlators at T = and 
the smeared correlators at T ~ 0.9Tc. For the lat- 
ter, the operators are spatially extended by Gaus- 
sian functions (/'(a?) oc exp (ajafp) with a — 0.16 
and 0.037 which are referred as "smeared" and 
"half-smeared" , respectively. 

First we present a result for the point opera- 
tors. Using the correlators at T = we test the 
applicability to finite temperatures by varying the 
number of degrees of freedom. Figure^shows the 
tmax dependence of the MEM results, where the 
fit ranges are restricted to < = 1 ~ tmax- The 
result at tmax = 16 can not reproduce even the 
position of the lowest peak, while tmax — 16 cor- 
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Figure 2. The spectral function for the smeared 
PS correlator at T ~ 0.9Tc. The default model is 
given by m(cj) = m]ji\iuj^. 

responds to the similar situation we encounter at 
T ^ 0.9Tc. This is one of the reasons why we 
give up to use the point operators to explore fi- 
nite temperatures. 

The MEM result with the smeared PS correla- 
tor at T ~ 0.9Tc are presented in Fig. 12 The low- 
est and next-lowest peaks are located at almost 
the same positions as at T = which is expected 
to have a similar result to that at T ~ O.QTc. This 
is encouraging for applications of the MEM to the 
systems at T > 0. Since the results are sensitive 
to the default model function (as input parame- 
ters), it is difficult to obtain quantitative results 
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Figure 3. Nterm dependences of mass and width 
by CCF for the smeared and half-smeared PS cor- 
relators. 

only with MEM. 

Based on the MEM results in Fig. [3 we assume 
multi-peak functions with the Breit-Wigner type 
forms, 
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as fitting functions and estimate the prior knowl- 
edge of the CCF. In Fig. El the fit resuhs for 
masses and widths of the lowest and next-lowest 
peaks are displayed, where the results for smeared 
and half-smeared correlators are plotted with dif- 
ferent symbols. Although with a sufficiently large 
Nterm (> 4) the results are stable, the differ- 
ences between the values for the smeared and 
half-smeared correlators are larger than the sta- 



tistical errors. Furthermore the results are rather 
sensitive to the prior knowledge. These results in- 
dicate that we need further investigations of CCF 
to obtain quantitative results with keeping sys- 
tematic uncertainties under control at finite tem- 
perature. 

Following our strategy, we have to perform 
these reliability checks with correlators at T = 
and to examine the t„iax dependence of the sys- 
tematic uncertainties. We also need further ex- 
amination of the smearing operators, because the 
suppression of high frequency part of the correla- 
tors may not be sufficient with the present smear- 
ing operators. This might be a reason that CCF 
does not work stably at T ~ 0.9Tc. These prob- 
lems must be clarified for reliable analysis of the 
hadronic correlators at finite temperature, in par- 
ticular in the plasma phase. 

This study was performed on the SR8000 at 
KEK (High Energy Accelerator Research Organi- 
zation) and the SX-5 at RCNP (Research Center 
for Nuclear Physics). T.U and H.M arc supported 
by the JSPS Research Fellowships for Young Sci- 
entists. 
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